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Abstract 

In this paper we investigate the model of opinion dynamics with anticon- 
formity on a complete graph. We show that below some threshold value of 
anticonformal behavior spontaneous reorientations occur between two stable 
states. Dealing with a complete graph allows us also for analytical treat- 
ment. We show that opinion dynamics can be understood as a movement of 
a public opinion in a symmetric bistable effective potential. We focus also 
on the spontaneous transitions between stable states and show that a typical 
waiting time can be observed. 

Keywords: Sznajd model, opinion dynamics, phase transition, stochastic 
resonance 



1. Introduction 

The word "revolution" came from the Latin "revolutio", which means 
"a turn around" - a fundamental change that takes place in relatively short 
period of time. Although, revolutions occurred many times in human his- 
tory, scientists still investigate what are possible reasons for a revolution. 
Particularly interesting empirical examples are recent protests in Egypt and 
Tunisia. In Egypt protests started on Tuesday, January 25 2011. These were 
the first protests on such a large scale to be seen in Egypt since the 1970s. 
What was the reason, why these protests begun on January 25? Sociologists 
claim that it was a domino effect - Egyptian revolution has been inspired by 
the successful revolution in Tunisia. Immediately the next question arise - 
what was the reason for the revolution in Tunisia? Again, sociologists and 
journalists claim that the answer is known - it has been sparked by the sui- 
cide of a young man who could not find a job and was barred from selling 
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fruit without a permit. There are of course many other examples of revo- 
lutions that took place in history. Usually, reasonable explanation is given 
by sociologists post factum. On the other hand, knowing the history one can 
always find some explanation. The problem is that we can hardly predict 
exact moment or place of a revolution. Why? Is it true that every overturn 
has a direct source? Physicist's experience shows that if a system supports 
two stable states, then fluctuations can cause random transitions between 
those two states [l], Recently such a spontaneous transitions, induced 
by a random noise, have been shown in a one-dimensional Sznajd model of 
opinion dynamics [3|. 

It should be noticed that the idea to introduce a random noise to the 
Sznajd model has been proposed already in the original paper from 2000 j^. 
However, in spite of the idea no systematic studies have been performed. 
In 2005 de la Lama et al. have included a social temperature effect in the 
model jsj. They have assumed that at every time step an agent followed 
the rules of the Sznajd model with probability p, while with a probability 
1 — p rules were not fulfilled (an agent adopted the opposite option than 
the one dictated by the rules) Such a modification has led to so called 
contrarians behavior j6|. It has been shown that for p > Pc the system was 
bistable with a probability density of magnetization having two maxima at 
m = m±. It has been claimed that for p > pc the system got ordered 
spontaneously selecting one of the stable solutions m = m±, while for p < Pc 
the system became monostable and disordered j^. Results obtained recently 
in (3| for one- dimensional system confirmed those obtained by de la Lama 
et al. Moreover, it has been shown that the finite system with random 
'contrarian-like' noise led to spontaneous transitions between two ordered 
stable states m = m± = ±1 [i*]. 

It should be recall here that the possibility of spontaneous transitions 
in the Sznajd model with noise has been also investigated by Wio et al. 
by including an external modulated field. This is well known that usually 
bistable systems response to a periodic external force in the presence of noise. 
This phenomena is known as a stochastic resonance and it has been widely 
studied in various systems over last 30 years [H, [sj. Wio et al. have shown 
that the system described by the Sznajd model with contrarians did not 
show conventional stochastic resonance, but so called system size stochastic 
resonance. 

In this paper we investigate different problem - spontaneous transitions 
not induced by any external force. We consider again generalized version of 
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the Sznajd model introduced in |3| in which two types of social response (con- 
formity and anticonformity) are possible. In ^ we have shown, performing 
Monte Carlo simulations for one-dimensional system, that for small values of 
anticonformity level, spontaneous reorientations occur between two opposite 
consensus states ('all spins up' or 'all spins down'). In this paper we consider 
the same model on a complete graph, which allows us for analytical treat- 
ment. Moreover, we show that results obtained for the complete graph are 
qualitatively different from those obtained for the one- dimensional system 
and more suitable to describe opinion dynamics in real social systems. 

We would like to clarify here that in econophysics and economy term 
'anticonformity' refers usually to contrarian behavior - a preference for tak- 
ing a position opposed to that of the majority. In we have used a term 
'anticonformity', that is usually used in social psychology jo], [lo[, to distin- 
guish between two commonly mixed up types of social response - 'anticonfor- 
mity' and 'nonconformity'. While anticonformity is similar to conformity in 
the sense that both (conformers and anticonformers) acknowledge the group 
norm (the conformers agree with the norm, the anticonformers disagree), 
nonconformity means independence. As long as we do not focus on partic- 
ular applications of the model (financial, social or any others) we can use 
interchangeably both terms 'contrarian' and 'anticonformists'. 

2. The model 

We consider a set of individuals, which are described by the binary 
variables S = 1 (t) or S* = — 1 (|) on a complete graph. At each elementary 
time step two individuals (denoted by f or J,) are chosen at random, and 
they influence a third randomly chosen individual (denoted by fl' or |l) in the 
following way: 

tt-l|~^ttlT with probability pi, no change with 1 — pi 
iifr^J-i-ll with probability pi, no change with 1 — pi 
ttff— T-ffJI with probability p2, no change with 1 — p2 
iiij-^iiil with probability p2, no change with 1 — p2. (1) 

First two processes correspond to conformity, and the next two describe anti- 
conformity. Above model is a particular case of the generalized Sznajd model 



introduced in It has been suggested in that significant parameter 



of the model is the ratio r = P2/P1 and changing absolute values of pi and 
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P2 influence only a time scale. Moreover, in the next section we show an- 
alytically that indeed stationary states depend only on the ratio r. In the 
previous paper [s^ we have decided to investigate the case in which p2 = 1 
and pi G (0, 1) is the only parameter of the model. In this paper we keep this 
assumption for Monte Carlo simulations, although we present also general 
analytic results for arbitrary values of pi,P2- 

In our previous paper [S], to investigate the model, we have provided 
Monte Carlo simulations with the random sequential updating mode and 
thus the time t has been measured in the Monte Carlo Sweeps (MCS) which 
consisted of N elementary updatings, i.e. At = As usually, we have 

defined a public opinion as a magnetization of the system: 



Monte Carlo results for the one-dimensional system have shown that for small 
values of anticonformity level p2 spontaneous reorientations occur between 
two extreme values of m = ±1. We have shown that there is a special value 
of anticonformity level p2 = p* below which the system stays for most of 
time in the one of two possible consensus states (m = ±1) and spontaneous 
reorientations occur randomly. Above the threshold, i.e. for p2 >P* system 
becomes monostable with magnetization fluctuating around zero [3|. 

In this paper we investigate the model on a complete graph, which allows 
us for analytical treatment. Moreover, we show that in case of a complete 
graph, spontaneous reorientations occur not between two extreme values of 
m = ±1, but between two stable solutions m = m± ioi P2/P1 = r < r* = 1/3. 
In other words, in the case of a complete graph we observe continuous phase 
transition at r*, while discontinuous transition in a case of one-dimensional 
lattice. 

In this paper we focus mostly on the spontaneous transitions between 
stable states, which were not investigated analytically up till now. Such a 
random transitions between two stable states are particularly important from 
social perspective to understand revolutions, protests and other fundamental 
social changes that take places in relatively short period of time. 

3. Analytical results for the infinite system 

In the case of a complete graph the state of the system is completely 
described by the magnetization m defined by Eq. ([2]) and takes the following 




(2) 



k=l 
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values: 

m = {-1,-1 + Ajv,-1 + 2A^,...,1}, (3) 

where = 2/N. 

Let us denote by A^^ number of spins 'up' and by number of spins 
down: 

N-^ + N^ = N 

N^-N^ = mN. (4) 

From above equations we can easily derive the formula for the probability of 
choosing randomly spin 'up' (concentration c of 'up' spins) as: 

We can now calculate the probability p+ that the concentration of 'up' spins 
increases c c+l/N {m ^ tji + Am) and the probability p_ that c decreases 
c — 7- c—l/N m — 7- m — Ajv). As seen from the definition of the model given by 
(II]) only 4 of 8 possible spin configurations are active tti; ttt; lit: ill- If cin 
inactive configuration is chosen, nothing happens. In the case of the infinite 
system, probabilities of choosing active configurations are the following: 

Pm = c^(i-c), 

Pm = c^ 

Pm = (1 - c)^c, 

PUi = (l-cf. (6) 
From (II]) transition probabilities can be calculated: 

Pm-^tn = PrPnh 
Ptn->m = P2-Pm^ 

PUt->U4 = Pi 'Put' 

(7) 

Probabilities p+ andp_ that the concentration of 'up' spins increases/decreases 
are equal: 



P+ = PrPm + P2-PUh 

P- = PrPm + PrPm^ (8) 
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and finally: 



P+ = (1 - cfP2 + C^(l - C)pi, 
P- = C^P2 + C(l - C)V- 



(9) 



After simple algebraic transformation we obtain the evolution equation: 



= c+^[~2{pi+p2)c^ + 3{pi+p2)c^-{pi + 3p2)c + p2], (10) 

where c' = c(t + A^) and c = c(t). We can easily calculate the fixed points 
c' = c of above equation: 



Cl = 
C2 = 

C3 = 

where useful quantity 



d 



Pi - 3p2 
P1+P2 

Pi - 3p2 " 
P1+P2 

Pi - 3p2 



Vd), 



(11) 



(12) 



Pi +P2 

For d < — T-r = P1/P2 > 1/3 the only real solution of equation flTU]) is 
c = Cl. For (i>0— !'r=pi/p2 < 1/3 we have two stable solutions c = C2 
and c = C3, while c = ci becomes unstable. From (ITT!) , using formula ([5]), we 
obtain immediately stationary values of magnetization in a very simple form 
(see Fig. [H): 



mst 



m± = ±Vd for r < 1/3 
mo = for r > 1/3. 



(13) 



It should be noticed that F = p^ — p_ can be treated as an effective force 
- p+ drives the system to the state 'spins up', while p_ to 'spins down': 



p_ = -2{pi + p2)c^ + 3(j9i + p2)c^ - {pi + 3p2)c + P2- (14) 
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Figure 1: Stationary values of magnetization m± as a function of the ratio r — p2/pi- 
Below the threshold vahie r = 1/3 there are two possible vahies of magnetization, above 
the threshold there is no majority in the system. From social perspective the system is in 
a stalemate state. 
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Figure 2: Effective potential y as a function of magnetization for pi — 1. As seen below 
threshold value p2 — 1/3 potential is bistable - there are two stable states at to = m± 



Therefore we can easily calculate the effective potential: 

y =^{Pl+ P2)C^ - {Pl + P2)C^ + ^{Pl + 3P2)C^ - P2C. (15) 

Using formula ([5]) we obtain the dependence between effective potential V 
and magnetization m. As seen from Figj2] below the threshold value r = 
P2/P1 < 1/3 potential is bistable - there are two stable states at m = m± 
given by equation ([T3|) . For r > 1/3 the only stable point is m = 0. From 
this perspective opinion dynamics can be understood as a movement of a 
public opinion m in a symmetric bistable potential V{m). Therefore we can 
expect random transitions between the neighboring potential wells, at least 
in a case of finite systems. 
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4. Analytical results for the finite system 



To investigate spontaneous transition between two metastable states, let 
us now derive the evolution equation for the finite system. Again, we start 
from calculating the probabilities p+ that the magnetization increases m — )■ 
m + Atv and p_ that magnetization decreases m — )• m — A^v- In the case 
of the finite system, probabilities of choosing active configurations are the 
following: 



Pm 
Pm 
Pm 

Pm - "jv ■ _ 1 ■ - 2 ■ ^"^^^ 

Again, using ([7]) transition probabilities can be calculated and finally: 

PiN^iN^ - l)N^+p2N^{N^ - 1)(A^ - 2) 
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A(A-l)(A-2) 
PiA;(A; - l)At + p2At(At - l)(At - 2) 



A(A-l)(A-2) 
Po = l-{p++p-). (17) 

Using formulas (jlj) we obtain: 

N, = l^N, iV, = i^iV, (18) 

and therefore we can rewrite transition probabilities p+,P- (flTI) as functions 
of m. 

Let P{m, t) denotes the probability density function of the magnetization 
m at time t. Time evolution of P{m^t) is described as usual by the master 
equation: 

P(m,t + At) =p+(m - A7v)P(m- Ajv,t) 
+p_(m + AN)Pim + A^, t) 
+ [1 - p+{m) - p_{m)]P{m,t), (19) 
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where = l/N, An = 2/N. 

Evolution equation allows to calculate numerically not only the asymp- 
totic states, but also the time evolution of the system, probabilities of transi- 
tions between two asymptotic stable states and probability density function 
(PDF) of waiting times (see the next section). 

For large, but finite systems (1 << < oo) we can solve the master 
equation analytically. Let us first rewrite it in the following form: 

P{m,t + At) - P{m,t) = 
= [p(m - A7v)P(m - Atv, t) + p(m + AAr)P(m + Aa^, t) - 2p{m)P{m, t)] 
-[F(m - A^)P(m - A^, t) - F{m + A^)P(m + A^, t)]/2, (20) 



where: 



. = ^ (21) 



and P = P+ — p_ is an effective force defined already in Eq. . 

For 1 << < oo above equation can be approximated by the following 
differential equation: 



which is a well known Fokker-Planck equation with diffusion coefficient: 

P+ +P- P1+P2 . 3p2 - Pi 2 /oqN 

P = ^^ = ^ + ^^- (23) 

and drift 



I 3P2-Pl . PI+P2 2 \ foA\ 

P+-P- = -[ ^ + — ^ — m ] m. (24) 



The stationary equation takes the following form: 



1 d NF 
which has a general solution: 

C f NF 
Pirn) = — exp / dm, (26) 

P J '^P 

where C is normalizing constant. Above equation has two qualitatively dif- 
ferent stationary solutions depending on the quantity d defined by ffT^ : 
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• For d = we obtain: 

P(m) = Cexp{-Nm^/A), (27) 

• For d> 0: 

P{m) = C(l - dm'y^l - dm^f'^~'-^^^'^ exp{Nm^ /2d). (28) 

P{m) has two maxima at mmax = = ±"\/^. It should be noticed here 
that rrimax is equal to rUst found in (fT3|) for the infinite system. 

We would like to stress here that above analytical formulas for the station- 
ary density probability function of magnetization P{m) have been derived 
under assumption that the system size is large i.e. iV >> 1. This allowed 
us to replace discrete by continues differential equation. However, obtained 
results agree very well with exact results, that can be obtained numerically 
from the evolution equation ( fT9|) . already for = 100 (see Fig. [3l). 

In the next section we compare our analytical and numerical results with 
Monte Carlo simulations and first of all investigate random transitions be- 
tween two stable states rUst = 'm±, that can be observed in a case of a finite 
system. 

5. Monte Carlo Simulations and comparison with analytical results 

We have performed computer simulations on a complete graph for p2 G 
(0, 1) for several system sizes - from = 50 to 200. As we have noticed in 
the previous section asymptotic behavior of our system is determined only 
by the ratio r = p2/pi between anticonformity and conformity. From social 
point of view the probability pi of conformal behavior is always greater then 
the probability p2 of anticonformal behavior, therefore we have decided to 
choose pi = I and p2 G [0, 1] (the same values of parameters pi,P2 has been 
chosen in ^). That means that for p2 = we have original Sznajd model. 

The time evolution of the public opinion, shown in Figure HI reminds a 
typical stochastic realization sampling the random waiting times in a sym- 
metric bistable potential, with stable states located at m = m± This is 
seen that values m± are p2-dependent and above certain threshold p2 > P2 
system becomes monostable. For p2 transitions between two stable state 
are rare, i.e. stability of states m = m± is high. With increasing aniconfor- 
mity level p2 transitions are more and more frequent - stability of m = m± 
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Figure 3: Stationary density probability function of magnetization from analytical formula 
([28|) (solid line) and obtained numerically directly from the master equation (|19|) (stars) 
for N — 25 (upper panel) and = 100 (bottom panel); pi — 1. 
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Figure 4: Time evolution of the public opinion for several values of anticonformity level p2 
and pi = 1 - for the low level of anticonformity spontaneous transition between two stable 
states occur. This behavior reminds a typical stochastic realization sampling the random 
waiting times in a symmetric bistable potential, with stable states located at m = m± 0] • 



decreases. This behavior is clearly seen from the dependence between PDF 
of public opinion and p2, which can be obtained from Monte Carlo (MC) 
simulations as well as analytical treatment presented in the previous section. 

In Figure [5]PDFs of public opinion for p2 = 0.2 for several lattice sizes are 
presented. As seen, numerical results obtained from ( fT9l) agree with Monte 
Carlo simulations, which is expected since we deal with a complete graph. In 
the next figure [6] we present probability density function of the public opinion 
P{m) depending on anticonformity level p2- We have shown only analytical 
results to make figure more legible but MC simulations give consistent results 
for the whole range of anticonformity level p2 G (0, 1). As seen from figured 
for small values of p2 the system is bistable with modes at m = m±. With 
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Figure 5: Probability density function of the public opinion P{m) divided by A^r = 2/N 
for several lattice sizes and anticonformity level p2 — 0.2. Numerical results obtained from 
equation agree with Monte Carlo (MC) simulations. 
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Figure 6: Probability density function of the public opinion P[m) for pi = 1 calculated 
numerically from Eq. ([T9| . 
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increasing p2 the distance and the hight of the separation between modes 
decreases. For certain threshold value p2 = P2 both modes approach m = 
and PDF of public opinion P{m) becomes unimodal. If we choose as an 
order parameter asymptotic absolute value of public opinion, we see that 
the system exhibits continuous phase transition at p2 = P2 = 1/3 (see Fig. 
[7]), as predicted analytically. Below the critical value of anticonformity two 
symmetric stable states exist. Due to fluctuations, spontaneous transitions 
between these two states are possible. Above the transition point p* there 
is no majority in the system, i.e. m fluctuates around zero. From the social 
point of view, community is in a stalemate state. It should be recall here 
that in the case of one-dimensional system we have also observed spontaneous 
transitions between two stable states. However, on contrary to the case of a 
complete graph, stable states were equal m = ±1 for any value of p2 < P2- 
Results on a complete graph are more reasonable for social systems - majority 
instead of unanimity is present for low level of anticonformity. Moreover, the 
level of the majority depends on p2. Similar results regarding the type of 
the phase transition and it's dependence on the topology were previously 
obtained by de la Lama et al., but spontaneous transitions between modes 
below the critical point were not considered in [Sj]. 

Let us now focus on those transition and measure the waiting time r that 
system spends in one of two stable states before jumping to the second one. 
We expect that waiting time r increases with the system size and decreases 
with P2. From Eq. [19]we are able to calculate the probability density function 
of waiting times. It occurs that the typical, i.e most probable waiting time 
exists (see Fig. |8]). As expected the typical waiting time increases with the 
system size and decreases with p2. For low values of anticonformity level 
P2 the most probable transition time Tmax grows almost exponentially with 
the system size N, which indicates that in the infinite system spontaneous 
transitions never occur (see Fig. |9]). 

6. Summary 

In this paper we have investigated a simple model of opinion dynamics 
with contrarians on a complete graph. Using Monte Carlo simulations and 
analytical calculations we have shown that below the critical level of anticon- 
formity the system is bistable spending most of the time in one two possible 
stable states. Above the critical value system becomes monostable and no 
majority exists in a system - so called stalemate situation. On contrary to 



16 



0.8 



0.6 



0.4 



0.2 



VI 1 

\i * 

X* 


o MC: N=50 
* MC:N=100 
V MC:N=200_ 
analytical 


o\| o 
oV 




a\ o 

\* 




^ 

\ ° 




\v 








1 1 i;^ 


owoooooo 







0.1 



0.2 



0.3 



0.4 



0.5 



Figure 7: Asymptotic absolute value of a public opinion as a function of anticonformity 
level p2- Continuous phase transition is clearly visible for p2 = P* = 1/3. Below this value 
there are two possible stable states. Due to fluctuations, spontaneous transitions between 
these two states are possible. Above the transition point p* there is no majority in the 
system. From the social point of view, community is in a stalemate state. Analytical 
result obtained for the infinite system agrees with Monte Carlo simulations. 
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Figure 9: Dependence between the system size N and the most probable waiting time 
T-max is shown in the left panel. In the right panel - dependence between the system size 
N and the probability P{Tmax)- As expected the typical waiting time increases with the 
system size and decreases with p2- 
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one-dimensional case the phase transition is continuous. From social perspec- 
tive, especially interesting results regard the spontaneous transitions between 
stable/metastable states. As shown here using analytical treatment, in a case 
of finite system spontaneous transitions between states can be observed and 
the typical waiting time exists. 

Random transitions between the neighboring potential wells caused by 
fluctuational forces were observed in various type of systems [2|. However, 
most of the papers concentrate on so called stochastic resonance - the re- 
sponse of s bistable system to a periodic external force in the presence of 
noise [H, Is], 0]. Yet in many nonlinear systems coherent transitions are not 
stimulated by an external force 12| and this is also probably the case of 
social systems. 
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